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Abstract. We present the results of a mid-infrared (5-16. 5//m) imaging spectroscopy survey of Young Stellar Objects (YSOs) 
and their surrounding environment in four low-mass star formation regions: RCrA, p Ophiuchi, Serpens and Chamaeleon I. 
This survey was performed using ISOCAM and its Circular Variable Filters (CVF) and observed 42 YSO candidates: we were 
able to obtain complete 5-16. 5yum spectra for 40 of these with a spectral resolving power of A//S.A ^ 40. A number of spectral 
. features were measured, most notably the 9.1pm silicate feature, the bending modes of both water and CO2 ices at 6.0 and 

15.2yum respectively and the well-known unidentified feature at 6.8yum. The strength of the unidentified feature was observed 
£ — , to correlate very well with that of the water ice bending mode and far less strongly with the CO2 ice bending mode. This 

suggests, in a manner consistent with previous observations, that the carrier of the unidentified feature is a strongly polar ice. 
Absorption profiles of the bending mode of CO2 ice are observed to show a significant long wavelength wing, which suggests 
that a significant fraction of the CO2 ice observed exists in a polar (H20-rich) phase. The sources observed in RCrA, p Oph and 
Serpens show similar spectral characteristics, whilst the sources observed in Cha I are somewhat anomalous, predominantly 
showing silicate emission and little or no absorption due to volatile ices. However this is consistent with previous studies of 
this region of the Cha I cloud, which contains an unusual cluster of YSOs. From comparisons of the strengths of the water ice 
CIh' an d silicate bands we detect an apparent under- abundance of water ice towards the sources in p Oph, relative to both RCrA 

and Serpens. This may be indicative of differences in chemical composition between the different clouds, or may be due to 
evaporation. Finally the C02:H20 ice ratios observed towards the sources in p Oph show significantly greater scatter than in 
^ ' the other regions, possibly due to varying local conditions around the YSOs in p Oph. 
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5~ 1 1- Introduction state transitions occur in the mid-IR, and a number of different 

molecules have been detected towards YSOs (see the review 
The circumstellar environment of young stars is difficult to by van Dishoeck & Hogerheiide fT999b . Spectroscopy of YSOs 
study, as in its early phases a young accreting and evolving has also connrm ed the presence of a number of molecul ar ices, 
star can be surrounded by sufficient gas and dust to cause a such as ices of watel . C o 2 a nd CH 3 OH (eg. Gibb et al. BOOTH 
large extinction. However, at visible and near-IR wavelengths B oogert et al. EOOOV In addition we also see both absorption 
this continuum extinction is smooth and contains little infor- and emission from refractory solids such as silicates and PAHs, 
mation on the nature of the surrounding gas and dust. The re- and infrared continuum emission. This continuum emission re- 
radiated energy is emitted mainly at wavelengths longer than sults from the reproC essing of stellar radiation by the dense dust 
2pm. Molecular vibrational-rotational gas-phase and solid- envelopes of the YSO, and this infrared excess is often used as 

a criterion for the identification of YSOs in photometric sur- 

Send offprint requests to: Richard Alexander, e-mail: veys. 

rda@ast.cam.ac.uk TT „ , , . 

* „ , , .. T „„ _„ . . . ... . . . Unfortunately, studies from the ground are constrained to 

Based on observations with ISO, an ES A project with instruments J 

funded by ESA Member States (especially the PI countries: France, work within atmospheric windows. In some cases, such as 

Germany, the Netherlands and the United Kingdom) and with the par- C °2, telluric absorption by the molecule at 15.2pm makes its 

ticipation of IS AS and NASA. study in astronomical objects impossible. Furthermore, the sen- 

** Present address: Institute of Astronomy, Madingley Road, sitivity of observations in the mid-IR is greatly reduced be- 

Cambridge, CB3 0HA, UK. cause telescope and atmospheric thermal emission greatly in- 
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crease the photon shot noise. Consequently, ground-based stud- 
ies longward of K band are generally restricted to observing 
relatively bright objects. The only solution to these problems is 
to move to an airborne or space platform, where even a modest 
sized telescope can outperform ground-based systems. 

Mid-IR spectra of young stars at high resolution were suc- 
cessfully taken with the Infrared Space Observatory (ISO, see 
Kessler et al. 1996 for an overview of the ISO mission) by var- 
ious investigators (eg. W33A by Gibb et al. pOOO; Elias 29 
by Boogert et al. 2000 1. In an attempt to push mid-IR ob- 
servations to even higher levels of sensitivity we made use 
of ISOCAM (Cesarsky et al. [T9961 and the Circular Variable 
Filters (CVF) to image 42 YSO candidates in four well-known 
low-mass star formation regions: RCrA, p Ophiuchi, Serpens 
and Chamaeleon I. The high sensitivity and good sample size 
was achieved by using the low CVF resolution and by observ- 
ing dense clusters of YSOs, so that ISOCAM's imaging capa- 
bility allowed simultaneous spectra to be obtained. It should 
be noted that the low spectral resolving power of the CVF 
(A/AA 40) essentially prohibits the detection of gas-phase 
spectral features; this is a study of the molecular ices and dust 
that lie around these young stars. However it is possible to 
draw a number of conclusions about the physics and chemistry 
around these YSOs, with information about their composition, 
thermal behaviour and evolution readily obtainable from such 
spectroscopic observations. 

2. Observations and Data Reduction 

2.1. Observations 

5-16. 5/mi imaging spectroscopy was performed on six tar- 
get fields in four low-mass star formation regions using the 
ISOCAM-CVF, on various dates between April 1996 and 
February 1997. This instrument made use of 2 variable width 
filters to take images at a number of different wavelengths, us- 
ing a 32x32 pixel gallium-doped silicon array. Exposure times 
of 2.1s per CVF step were used. Each was exposed repeatedly 
and the images were sequenced (in time) in decreasing order 
of wavelength, resulting in a total exposure time of ~2500s per 
target region for each CVF set. This procedure produced 2 dat- 
acubes for each target region. The first covered the wavelength 
range 5.0-9.4/mi (CVF1) and the other covered 9.3-16.3//m 
(CVF2), with a mean spectral resolving power of A/AA =* 40. 
A brief summary of the characteristics of the observed regions 
is given below. 

2.1.1. RCrA 

The RCrA cloud is the name given to a molecular cloud of mass 
~ 120Mq (Harju et al. 1993 1 around the variable star R Coronae 
Australis, at a distance of 130pc (Marraco & Rydgren 198T}- It 
is a region of low-mass star formation located approximately 
18° below the galactic plane (Olofsson et al. [19991 , where 
it suffers little foreground obscuration. The previous K-band 
studies of Taylor & Storey JT9841 and Wilking et al. ([1997) and 
the sub-millimetre survey of Harju et al. ( 1993 ) revealed a num- 
ber of embedded objects and other objects classified as YSOs. 



One field of view was observed in the RCrA cloud, using the 
6" pixel field of view (pfov). It was centred on 19 h 01 m 45 s .l, 
-36°58'34" 1 , immediately to the south-west of RCrA itself. 

2.1.2. p Ophiuchi 

The p Ophiuchi dark cloud complex is a complicated struc- 
ture of several large molecular clouds near the bright star p 
Ophiuchi. Its relatively high galactic latitude (b-\l°) com- 
bined with its relative proximity (t/=160pc, Whittet II 974> re- 
sult in comparatively little foreground obscuration, making it 
an ideal region in which to study star formation. Initially iden- 
tified by Grasdalen et al. d!973> . who identified 41 sources in a 
K-band map, it has since been studied extensively both in the 
infrared (Elias fT978l Greene & Young[T332| Strom et al. [T995l 
Abergel et al. 1996, Bontemps et al. 2001 1 and at millimetre 
wavelengths (Loren et al. 119901 Andre & Montmerle I1994> . 
In particular, Loren et al. (1990 1 showed that the cloud con- 
sists of a number of smaller sub-structures. They identified 6 
dense cores (labelled A-F) and observed significantly higher 
clustering of YSOs in three of these cores (A, B, and E/F). Two 
fields were observed in p Oph, both with the 6" pfov: one cen- 
tred on 16 h 26 m 21 s 4, -24°23'59" (field p Oph A) and one at 
16 h 27 m 23 s :7, -24°40'39" (field p Oph E). The positions corre- 
spond (approximately) to cores A and E/F above. 

2.1.3. Serpens 

The Serpens molecular cloud was initially identified by Strom 
et al. ( 1976 1. It has since been very well surveyed, both in the 
infrared (Churchwell & Koornneef 1986 Gomez de Castro et 
al. fT988l Eiroa & Casali IT9921 Giovannetti et al. fT998l Kaas 
1999 1 and the sub-millimetre (Casali et al. 1993 1, with as many 
as 163 sources identified in the 6'x5' K-band map of Eiroa 
& Casali ( 1992 1. The cloud itself extends over a region of 
>20'x >15' in CO line emission (Loren et al. 1979 1 and is at a 
distance of ^260pc (Straizys et al. H996i . In addition to survey 
imaging, many individual objects in the cloud have been stud- 
ied in detail: Eiroa & Casali ( 1989 1 studied the multiple outflow 
source SVS4; SVS2 has been shown to illuminate the Serpens 
Reflection Nebula (Gomez de Castro et al. 1988 Huard et 
al.Q997i; SVS20 has been shown to be a close binary system 
(Huard et al. [T997l . Two fields were observed in Serpens. Field 
Ser B, taken with the 6" pfov, was of the region around SVS2 
and SVS20 and was centred on 18 h 29 m 57.6, 1°12'42". Field 
Ser A was a higher resolution field (3" pfov) of the multiple 
outflow source SVS4, centred on \&29 m 5T2, 1°14'25". 

2.1 .4. Chamaeleon I 

The Chamaeleon dark cloud complex is a complicated struc- 
ture consisting of 3 large molecular clouds (designated Cha 
I,II,III by Hoffmeister 1963 1 and a number of smaller clumps 
and globules. Again, the relative proximity (ii=160pc, Whittet 
et al. 119971 and high galactic latitude (b ~ -16°) of the clouds 
result in little foreground obscuration and make them ideal for 

1 Please note that all coordinates given are J2000.0. 
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star formation study. The region has been well-surveyed at mil- 
limetre wavelengths (Mattila et al. 1989 Henning at al. 119931 
and more recently in the near- (Kenyon & Gomez 2001) and 
mid-infrared (Persi et al. 2000). A single field was observed 
in the Cha I cloud, with the 6" pfov, centred on ll h 09 m 39 s :4, 
-76°35'2". The corresponds (approximately) to a previously 
observed dense core lying to the north of HD97300 (Jones et 
al. 113551 Persi et al. rT333l . 

2.2. Data Reduction 

Initial data reduction was performed using the auto-analysis 
package from the ISO data pipeline (OLP VI 0.0, Ott et 
31.03)97, Blommaert et al. 2001). This corrected the images 
for changes in the spectral response of the detector, as well as 
performing flatfielding and flux calibration. It should be noted 
that column 24 on the array was dead: this was corrected for by 
means of a simple linear interpolation from the adjacent pixels. 
Sources were identified "by eye" and everything representing 
a local maximum in intensity in more than 75% of the frames 
was treated as a source. For each source identified the flux in 
each frame was evaluated as the total flux contained within a 
18"xl8" (3x3 pixel) aperture 2 less the sky flux (evaluated as 
the mean flux per pixel over a designated sky region). This pro- 
duced a spectrum for each source from both filter sets. These 
spectra were smoothed to match the resolution of the CVF and 
the 2 "halves" taken from the 2 filter sets were then joined to- 
gether. Unfortunately it was necessary to discard the data from 
positions 70-80 on CVF1 (8. 94-9. 43pm), as the array "mem- 
ory" (a noted problem with ISO, see Blommaert et al. 2001 ) 
rendered these images useless. Consequently the spectra con- 
tain a small gap in their wavelength coverage at this point, from 
8. 94-9. 33pm. The spectra were also corrected for a systematic 
under-estimation of flux at long wavelength, which arose from 
the use of a fixed-size aperture to measure the total flux from 
diffraction-limited images. 

2.3. Spectral Features 

An example of a typical spectrum is shown in Fig[J The fol- 
lowing spectral features were identified in the spectra: 

- An absorption feature at 6pm attributed to the bending 
mode of H2O ice (Gerakines et al. 1995 ). 

- An absorption feature at 6.8pm, the well-known unidenti- 
fied feature (Schutte fT337l 

- A weak absorption feature at ~7.6pm which was not re- 
solved at the low spectral resolution of the CVF. It is ten- 
tatively identified as the bending mode of CH4 at 7.67pm 
(Hudgins et al. 119931 Boogert et al. 119961 Boogert et 
al. 1998 ) but due to the low resolution of the CVF and the 
weak nature of the observed feature this identification is 
somewhat uncertain. 

- A broad feature centred at ~9.7pm, seen in both emis- 
sion and absorption, identified as the Si-O stretching mode 
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Fig. 1. A typical spectrum, with prominent spectral features an- 
notated. 



of "astronomical silicate" (Draine & Lee 1984, Simpson 
[T33T1 O'Donneli rrffHl . 

- A broad absorption feature at ~13/im identified as the li- 
bration (or hindered rotation) mode of H2O ice (Gerakines 
at al. fT3351 . 

- An absorption feature at 15.2pm attributed to the bending 
mode of CO2 ice (Ehrenfreund et al. 119961 . 

Other spectral features were observed in individual sources, 
such as the apparent 10.5pm emission feature in RCrA IRS5 
(see Fig0 or the apparent 14pm absorption feature in HH100- 
IR (see Fig0J. These rarer features have been noted in the re- 
sults tables where they occur, but the goal of this study was to 
look for general trends rather than individual curiosities. With 
this in mind, and also considering the large number of spec- 
tra obtained, it was decided to fit only the strongest and most 
commonly occuring features highlighted above. 

2.4. Fitting the Spectra 

As many of these spectral features overlap, independent de- 
terminations of their strengths were not possible, so all of the 
features in the spectra were fitted simultaneously. This was 
achieved by fitting a spectrum of the form: 



F y (A) = C(A) x exp 



• 2^ Siai(A) 



(1) 



2 A 15"xl5" (5x5 pixel) aperture was used for the 3" pfov frames 
offieldSerA. 



where C(A) is the continuum, s,- represents the absorbing col- 
umn density of each species i and ai(A) are the absorption co- 
efficients for each species i. This essentially treats the absorb- 
ing material as a cold absorbing screen. In order to account for 
possible silicate emission, the depth of the silicate feature, ssu 
was allowed to run negative. Note however that this assumes 
that the same profile applies in both absorption and emission, 
so could result in errors if both are present. The absorption 
coefficients for H2O ice and CO2 ice were taken from the 
Leiden Observatory Database of Interstellar Ices (Gerakines et 



4 



R.D. Alexander et al.: ISOCAM-CVF Spectroscopy of YSOs 



al. 119951 Ehrenfreund et al. 1996 1 and those coefficients eval- 
uated at 10K were used, as thermal variations are negligible at 
the spectral resolution of the CVF. It should be noted, however, 
that the depths of the 6fim and 13yum water ice bands did not 
correlate as expected, so they were fitted independently. The 
unidentified feature at 6.8yum was fitted using an asymmetric 
Gaussian profile: a Gaussian profile with two different widths 
shortward and longward of the peak. The best-fitting profile for 
each 6.8jt/m feature was measured, and the widths were then 
fixed as the means of the widths measured in all of the ob- 
served profiles (FWHM=500nm). The poorly resolved 7.6/um 
feature was fitted using a Gaussian profile with a width equal 
to that of the CVF resolution element at that wavelength. A 
number of profiles exist in the literature for the silicate feature, 
obtained both theoretically and observationally, so the profile 
which provided the best fit to the observed data was used: 
profile 2 from Simpson ( 1991 1. However the profiles obtained 
by Simpson d!991> from IRAS observations extend shortward 
only to ^7.5/im, and so shortward of this wavelength the op- 
tical constants for "circumstellar silicates" from Ossenkopf et 
al. (I1992t were used. It should be noted that the profiles from 
Simpson (I1991> . Ossenkopf et al. ( 1992 1 and Draine ([1985> are 
almost identical away from the 9.7/vm peak: only the shape of 
the absorption band is dependent on the choice of optical con- 
stants. 

The fitting procedure was iterative, with repeated fitting 
converging to a final profile. The fitting algorithm took the 
form: 

1. Fit initial continuum and divide out. The continuum was 
fitted as a 2nd-order polynomial through points at 5.6, 7.9 
and 16.1/mi, so chosen as they best constrained the result- 
ing fits. 

2. Measure the strength of each spectral feature at a number 
of points around the peak. 

3. Use these strengths and the optical constants for the differ- 
ent species to divide out the features from the spectra. 

4. Fit a new continuum to this "featureless" spectrum. 

5. Iterate this procedure until the fit converged: typically 10- 
15 iterations were sufficient. 

It should be noted that this allowed the continuum to drift away 
from the observed spectrum in some cases, due to absorption 
across the entire observed wavelength range. This is because 
the normalised silicate opacity over the range 5-7.5^im is ap- 
proximately constant at ^0.15, and towards the long wave- 
length end of the spectra the 18/mi bending feature begins to 
"cut in". (The normalised opacity at 16/im is approximately 
0.3.) The effect of this can be seen in Fig|2 although the con- 
tinua are seen to deviate significantly from the observed spec- 
tra, appearing to diverge from the observed spectra in some 
cases, this is primarily due to silicate absorption across the en- 
tire wavelength range. With only the silicate features included, 
these "continua plus silicate" fits are seen to follow the ob- 
served spectra well. As a result, any possible systematic errors 
arising from the manner of the continuum fit will only affect the 
evaluation of the silicate optical depth; the measured strengths 
of the narrower features are not affected by the choice of con- 
tinuum. 
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Fig. 2. Examples of successful fits to GY265 in p Oph (upper) 
and RCrA IRS5 (lower): fitted spectra are shown as solid lines 
and continua as dashed lines. The dotted lines are the "contin- 
uum plus silicate" fits given by C(A) x exp (-ssiffsi(^))- 

On inspection of the residuals from this fitting procedure 
a further feature was observed: a broad feature centred on 
^11.2/im, seen in both emission and absorption. This has a 
number of possible identifications, such as crystalline silicates 
(Bregman et al. nWfl Campins & Rvan [T9^9l . the Unidentified 
(PAH) Band at 11.3/im (Boulanger et al. 1996 1, or a shoulder 
on the silicate profile. Given the variation in possible identifica- 
tions, it was decided to fit this feature independently using an 
asymmetric Gaussian profile, once again taking the widths to 
be the mean of those observed (FWHM= 1700nm). This feature 
was then included in the algorithm and the fitting procedure re- 
peated: examples of successful fits are shown in Fig|2] Another 
similar residual around 8.5/mi was also observed, but the prop- 
erties and strength of this were greatly affected by the choice of 
continuum point around 8/im. Whereas the llyum feature was 
robust against small changes in the continuum this apparent 
8.5^m feature was not, so it was not included in the fitting pro- 
cedure. 

The 15.2/im CO2 ice feature is an interesting one and has 
been used in the past as a diagnostic of the ice environment. 
Gerakines at al. ([T999 1 found that two distinct phases of CO2 
ice (polar and non-polar) show significantly different absorp- 
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tion profiles. The non-polar phase is characterised by a double- 
peaked profile, arising from the two-fold degeneracy of the 
bending mode, whereas the polar (I^O-rich) phase shows a 
significant long wavelength wing. The spectral resolution of 
the CVF is too low to resolve the double-peak structure, but 
the CO2 profiles were observed to vary noticeably from source 
to source. As a result it was decided to let the CO2 profile vary, 
so it was fitted using an asymmetric Gaussian profile with the 
widths as variable parameters. 

There are two forms of error associated with this fitting 
procedure. Firstly there are the ordinary statistical errors as- 
sociated with the fits, evaluated by combining the statistical er- 
rors on the data points used to fit the features and the intrinsic 
uncertainties in the individual feature fits. These are typically 
around 5-10% for the stronger features (silicate, H2O, CO2) 
and around 15-20% for the less well constrained features, al- 
though in the case of some poor S/N spectra they are much 
larger. In addition to this, there are systematic errors associated 
with poor fits. These are usually due to the presence of "extra" 
features in the spectra, or because of poorly fitting feature pro- 
files, and are less easy to quantify. Consequently, the errors re- 
ported in the results tables are the statistical errors only; where 
fits are flagged as poor additional systematic errors also apply. 



2.5. Dealing with poor fits 

Generally the procedure described above proved successful, 
but in some cases it led to poor fitting of certain features. 
This was primarily due to the manner in which the continuum 
was fitted, as shown in the examples in Fig|3] The continuum 
was fitted as a 2nd-order polynomial through three sections 
on the spectrum expected to show the least absorption, as de- 
scribed above. However, constraining the continuum through 
these three points led in some cases to poor fitting elsewhere, 
particularly at the short end of the CO2 ice band and immedi- 
ately shortward of the silicate feature. 

In such cases, the narrower features at the short end of the 
spectrum (6.0, 6.8 and 7.6^/m) were fitted in exactly the same 
way as before, only this time using a linear continuum through 
the points at 5.6 and 7.9/vm. As absorption due to the silicates 
is approximately constant over this wavelength range this re- 
moved complications due to the quadratic continuum fit. When 
tested on spectra which were fitted well, such as those in Fig|2] 
this produced results which were the same as before, within the 
fitting errors. Also, the manner in which the continuum was fit- 
ted often led to the CO2 feature fitting poorly, as seen in Fig|3 
As a result it was decided to measure the CO2 absorption in all 
the sources by fitting a linear continuum from 14.5-16. lyum, as 
the silicate absorption does not vary much over this range and 
the CO2 is not blended with any other spectral feature. This 
procedure was only applied to these narrower features, how- 
ever, as the broad silicate, 1 1 and 13/mi bands were deemed to 
be too heavily blended to deconvolve them meaningfully over 
a short section of the spectra. In addition, all features fitted in 
this manner are flagged as such in the results tables. 
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Fig. 3. Examples of poor fits (solid lines) from 5-8/vm and 
around the CO2 ice feature. 

2.6. Spectral Index 

In order to take a standard measurement of the SED, the spec- 
tral index a was evaluated as a traditional IR spectral index 
(eg. Kenyon & Hartmann 1995i : 

log A 2 F, h - log A t F Al 

a = i ] i i W 

log A 2 - log A x 

This sets a — as the index of a flat vF v and a — -3 as the in- 
dex of a Rayleigh-Jeans spectrum. Equivalently, class I YSOs 
have a > 0, class II have -3 < a < and class III a ~ -3, 
using the conventional definition of spectral class (Lada 1987 1, 
although more recent studies (eg. Andre & Montmerle 1994; 
Greene et al. 1994) have found that the class II/III border lies 
at a higher value {a —1.5). Where previous K-band data ex- 
ist in the literature the wavelengths were taken as A 2 = 14.0/im 
(our data) and Ai = 2.2/im (K-band data from references in 
Table [0 and this was applied to both the observed spectrum 
and the calculated continuum for each source. (Where no K- 
band data exist, for sources not previously observed or with 
uncertain identifications, A\ = 8.0^m was used.) This second 
("continuum") measurement of the spectral index is important, 
as many of the sources suffer significant absorption (includ- 
ing foreground absorption) and a straight application to the ob- 
served spectrum could result in sources being mis-classified in 
cases where deep absorption is present. In cases where the flux 
suffers little extinction the two give the same results, but where 
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deep absorption is present the continuum spectral index derived 
is greater than that obtained from the observed spectrum. In 
such cases the continuum spectral index provides a truer mea- 
surement of the shape of the underlying continuum. 

3. Results 

3.1. Detected Sources 

A total of 43 sources were detected: a complete list is presented 
in Table [0 The 6" pfov used here results in coordinates which 
are accurate only to ±3—4" at best, so where sources have been 
identified with previous studies the coordinates presented are 
those from previous observations, where higher spatial resolu- 
tions or better S/N have resulted in more accurate positions. 
The coordinates evaluated here are only presented for sources 
which represent new or uncertain detections. In the 3 cases 
where sources have been identified as unresolved doubles the 
coordinates presented are those of the brighter source. Of the 
43 sources detected 6 were not identified with detections from 
previous surveys (see Sect. B.l.H and 2 (Ser A 7 and Ser B 1 1) 
were determined to be different images of the same source, as 
the two Serpens fields overlap slightly. Consequently we were 
able to obtain complete 5.0-16.3/mi spectra for a total of 40 
sources, along with a further 2 partial spectra 3 : a complete at- 
las of these 42 spectra is presented in AppendixlAl 

3.1 .1 . Ghost Images and New Detections 

As noted by Blommaert et al. (2001 1 ghost images and stray- 
light are a significant problem when using the CVF, so care 
must be taken when identifying sources. Ghost images arise in 
two ways and are found either around the true image or directly 
opposite it (relative to the optical axis). Potential ghost images 
were identified by looking at the spectra of the objects, their 
flux densities and their positions relative to the instrument op- 
tics. In this manner 7 objects initially identified as sources were 
reclassified as ghost images (4 in p Oph E, 2 in Cha I and 1 in 
RCrA). This leaves 6 "new" detections in our survey, which we 
now address: 

RCrA: 3 sources in RCrA are identified as new detections: 
sources RCrA 2, 6 and 7. Source RCrA 2 lies close to RCrA 
1, separated by ~30", but shows sufficiently different spec- 
tral characteristics that it cannot be a ghost of this source. It 
is therefore considered to be a real source, probably a YSO. 
Both RCrA 6 and 7 lie in a region of nebulosity between the 
bright sources HH100-IR (RCrA 5) and RCrA (which lies 
just outside the image frame). These are not considered to 
be ghost images, but may well represent peaks in the nebu- 
lous emission rather than true YSOs. It should also be noted 
that this region of the RCrA cloud was not included in the 
ISOCAM survey of Olofsson et al. i 19991 1. due to detector 
saturation, so it is not unreasonable to expect new detec- 
tions here. 

3 These partial spectra were the result of various observational prob- 
lems, such as sources lying near to the edge of the frame or lying on 
the bad column. 



p Oph E: A large number of ghost images were found here 
and after they were removed 1 unidentified source re- 
mained: p Oph E 4. However, this is the faintest source 
in the survey, with a flux rising to just 0.17Jy at 16^/m, so 
is probably a peak in the nebulosity rather than an "new" 
YSO. 

Ser B: 2 sources in Ser B were not identified with previous 
observations: Ser B 4 and 7. Both these sources lie near to 
SVS20 and so could be ring-like ghost images associated 
with it. However they also lie in a previously identified re- 
gion of nebulous emission (Kaas 1999 1 and so may instead 
be peaks in this nebulosity. It is impossible to resolve this 
ambiguity completely, but in either case these 2 sources are 
probably not "new" YSOs. 

3.1.2. Background Stars 

Of the 36 sources identified with previous studies 3 have pre- 
viously been identified as stars which are not associated with 
the star-forming clouds: CK2 (Ser B 12) by Chiar et al. JT9941 
and Casali & Eiroa Jl993b : Ser B 1 and Ser B 1 1 / Ser A 7 by 
Giovannetti et al. (I1998> . 

CK2 is the only one of these 3 sources which has been 
identified as a background field star (probably a background 
supergiant, Casali & Eiroa 1996 1. Unfortunately it fell directly 
on the bad column in the array here and so no useful spec- 
troscopic data regarding CK2 was obtained. Both of the other 
two sources in Serpens show extremely deep absorption fea- 
tures due to ices and silicates and both have very red SEDs. 
Taking the K-band magnitudes of Giovannetti et al. ( 1998 ) and 
assuming the SEDs of these stars to be Rayleigh-Jeans spec- 
tra (F v oc A~ 2 ) results in predicted mid-IR flux densities of 
~0.1mJy for both sources, approximately 1000 times less than 
what is observed. Consequently it seems unlikely that these are 
background sources, so we interpret them to be deeply embed- 
ded objects: with this interpretation these 2 sources are 2 of the 
3 most deeply embedded objects in this survey. 

3.2. Spectral Fitting 

The spectra of the 43 sources identified in Table ^ were fitted 
using the method described above: the results of this procedure 
are presented in Tableland the fits and continua obtained are 
presented in AppendixlAl The fitting procedure is both robust 
and unambiguous, converging to a good fit in most cases, but 
a few weaknesses exist. Firstly, as discussed in Sect. 12.31 the 
procedure only measures the 7 spectral features observed in al- 
most all of the sources and consequently does not fit some rarer 
features, such as an apparent 14yum absorption band observed 
in source p Oph A 2. Further, the profile fitted to each feature 
(except the CO2) is not allowed to vary from source to source. 
While this is broadly valid, a few exceptions led to some poor 
fits, such as the broadened 6.0 and 6.8/im bands in SVS2 (see 
FigEJ, or a handful of sources which show broadened silicate 
features. However, given the low spectral resolution of the CVF 
and the variation in the spectra over the large number of sources 
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Table 1. Detected sources: references indicate origin of specified coordinates. 



Source 


o^nnn (h m s) 


<5innn i ' ") 


Previous ID a 


Reference" 


Notes 


RCrA 1 


19 01 41.5 


-36 58 29 


IRS2, TS13.1 


2 




RCrA 2 b 


19 01 41.5 


-36 58 59 




17 




RCrA 3 


19 01 41.7 


-36 59 54 


Hff2 


1 




RCrA 4 


19 01 48.0 


-36 57 19 


IRS5, TS2.4 


2 




RCrA 5 


19 01 50.7 


-36 58 07 


HH 100-IR, IRS1,TS2.6 


2 




RCrA 6 b 


19 01 50.9 


-36 57 37 




17 


Very close to RCrA 


RCrA 7 b 


19 01 52.4 


-36 57 37 




17 


Very close to RCrA 


p Oph A 1 


16 26 17.3 


-24 23 49 


SKS9, ISO-Oph 21 


6 




p Oph A 2 


16 26 18.8 


-24 24 21 


SKSll,ISO-Oph 26 ??? 


17 


Uncertain ID 


p Oph A 3 


16 26 21.5 


-24 23 07 


GSS30, GY6, SKS12, ISO-Oph 29 


6 




p Oph A 4 


16 26 23.7 


-24 24 40 


GY21, SKS16, ISO-Oph 37 


6 


Falls on bad column 


p Oph E 1 


16 27 21.6 


-24 41 43 


GY252, SKS36, IKT43, ISO-Oph 132 


6 




p Oph E 2 


16 27 24.8 


-24 41 02 


ISO-Oph 137 


6 




p Oph E 3 


16 27 26.9 


-24 39 23 


GY262, IKT53, ISO-Oph 140 


6 




p Oph E 4 b 


16 27 27.1 


-24 41 31 




17 


Very faint 


p Oph E 5 


16 27 27.2 


-24 40 49 


GY265, IKT54, ISO-Oph 141 


6 




p Oph E 6 


16 27 28.2 


-24 39 32 


GY269, IKT57, ISO-Oph 143 


6 




Ser A 1 


18 29 56.6 


+ 1 12 56 


SVS4/2, GEL5 


11 


Very faint 


Ser A 2 


18 29 56.7 


+ 1 12 35 


SVS4/3 


11 


Very faint 


Ser A3 


18 29 57.5 


+ 1 12 57 


SVS4/5 and 
SVS4/6, GEL9 


11 


Unresolved double 


Ser A 4 


18 29 57.9 


+ 1 12 25 


SVS4/7 


11 


Very faint 


Ser A 5 


18 29 57.9 


+ 1 12 34 


SVS4/8 


11 




Ser A 6 


18 29 57.9 


+ 1 12 43 


SVS4/9 and 
SVS4/10, GEL 11 


11 


Unresolved double 


Ser A T 


18 30 00.1 


+1 13 03 


GCNM130 


13 


Off frame 7.0-9.3yt/m 


SerB 1 


18 29 52.9 


+ 1 14 55 


GCNM53 


13 




SerB 2 


18 29 55.7 


+ 1 14 31 


CK9, GEL4, EC74, GCNM76 


13 




SerB 3 


18 29 56.8 


+ 1 14 46 


SVS2, CK3, GEL6, EC82, GCNM87 
and GEL8, EC86, GCNM93 


13 


Unresolved double 


Ser B 4 b 


18 29 57.2 


+ 1 13 28 




17 


Peak in nebulosity? 


SerB 5 


18 29 57.6 


+ 1 15 31 


GCNM100 


13 




SerB 6 


18 29 57.7 


+ 1 14 05 


SVS20, CK1, GEL 10, EC90, GCNM98 


13 


Double Object 


Ser B 7 b 


18 29 58.1 


+ 1 13 24 




17 


Peak in nebulosity? 


SerB 8 


18 29 58.2 


+ 1 15 21 


CK4, GEL12, EC97, GCNM106 


13 


Falls on bad column 


SerB 9 


18 29 58.7 


+ 1 14 26 


EC103, GCNM112 


13 




Ser B 10 


18 29 59.2 


+ 1 14 08 


CK8, GEL13,EC105, GCNM119 


13 




SerB ll c 


18 30 00.1 


+ 1 13 03 


GCNM130 


13 




Ser B 12 


18 30 00.5 


+ 1 15 20 


CK2, EC 118, GCNM136 


13 


Peaks on bad column 


SerB 13 


18 30 02.1 


+1 13 59 


CK7, EC125, GCNM154 


13 


Very faint 


Chal 1 


11 09 23.3 


-76 34 35 


Cl-6, ISO-Chal 189, KG82 


15 


Bad column in CVF1 


ChaI2 


11 09 29.2 


-76 33 30 


ISO-Chal 192, KG87 


15 


Out of frame in CVF1 


Chal 3 


1 1 09 42.6 


-76 35 01 


Cl-25, ISO-Chal 199, KG93 


15 




Chal 4 


11 09 46.9 


-76 34 49 


CI -24, ISO-Chal 204, KG97 


15 




Chal 5 


11 09 54.1 


-76 34 26 


Cl-5, ISO-Chal 223, KG109 


15 




Chal 6 


11 10 00.8 


-76 34 59 


WW Cha, CI -7, ISO-Chal 231, KG116 


15 





a Names and references from: 1) Marraco & Rydgren fiWTl (Ha); 2) Taylor & Storey (1984} (TS, IRS); 3) Grasdalen et al. fT973l (GSS); 
4) Greene & Youn g f!992l (GY); 5) Strom et al. 119951 (SKS ); 6) Bontemps et al. OOOTl (ISO- Oph); 7) Imanishi et al. i200Tl (I KT); 8) 
Strom et al. (19761 (SVS); 9) Churchwell & Koornneef (T9861 (CK); 10) Gomez de Castro et al. (T^l (GEL); 1 1) Eir oa & C asali (T5H5J 
(SVS4/?); 12) Eiroa & Casali (EC); 13) Giovannetti et al. fEMl (GCNM); 14) Hyland et al. (TW2) / Jones et al. (19851 (CI-?); 15) 

Persi et a l. 1 200 0) (ISO-Chal); 16) Kenyon & Gomez i200fl (KG); 17) this paper. 

b See Sect. l3.1.T1 for a discussion of the sources not previously identified. 

These two sources from the two different Serpens fields were determined to be different images of the same source. 



Table 2. Fitting results: errors are given in parentheses (ie. 7.84(16) «-> 7.84+0.16); source names correspond to those in Table[T] 





8/jm 


8yum 








Equivalent Line Widths (nm) 




Source 


Cont. 






Source 


flux 

(Jy) 


continuum 
flux (Jy) 


Tsi" 




6/jm 
H 2 


6.8pm 


7.6pm 
CH 4 ?? 


13/jm 
H z O 


15.2yum 

co 2 


a 


a 


Group b 


Notes 


RCrA 1 


7.84(16) 


7.79(16) 


-0.02(1) 


-0.06(1) 


178(8) c 


17(l) c 


17(3) c 


247(25) 


168(7) 


0.56 


0.58 


b 


Strong 1 l//m emission 


2 


0.95(2) 


1.24(3) 


0.72(3) 


-0.27(4) 


203(8) 


57(4) 


17(3) 


0(0) 


83(4) 


-0.96 f 


-1.35' 


a 


Silicate too wide at long 
A, probable 13/jm error 


3 


0.27(2) 


0.28(2) 


0.22(1) 


-0.16(2) 


147(15) c 


14(3) c 


0(0) c 


718(76) 


0(0) 


-1.11 


-0.96 


b 


Tsi probably too small, 
possible absorption bands 
at 13pm, 14.25/um 


4 


1.74(110) 


2.28(144) 


0.71(3) 


-0.07(1) 


271(10) 


81(6) 


6(1) 


586(59) 


253(11) 


1.23 


1.43 


a 


Possible 10.5/um emission 


5 


22.46(44) 


36.91(74) 


1.35(5) 


-0.23(3) 


206(8) 


67(4) 


14(2) 


1535(155) 


129(5) 


1.30 


1.79 


a 




6 


0.94(3) 


1.36(4) 


0.95(4) 


-0.15(2) 


204(7) 


61(4) 


11(2) 


73(7) 


97(4) 


-0.40 f 


-0.42 f 


a 


Under-estimates 13/um 


7 


1.20(27) 


1.62(37) 


0.74(3) 


-0.08(1) 


191(9) 


86(6) 


18(3) 


0(0) 


95(4) 


-0.1 8 f 


-0.37 f 


a 


Possible 10.5/um emission, 
under-estimates 13pm 


p Oph A 1 


0.42(2) 


0.85(4) 


2.11(11) 


-0.27(4) 


27(3) 


72(7) 


19(3) 


927(94) 


57(3) 


-0.70 


-0.20 


a 


Poor 5-8pm fit 


2 


0.21(2) 


0.54(6) 


2.68(19) 


-0.70(10) 


207(14) 


253(53) 


79(9) 


2611(319) 


0(0) 


0.34 


1.36 


a 


Definite 14pm band 


3 


6.03(7) 


8.99(10) 


1.05(4) 


-0.20(3) 


107(6) c 


0(0) c ?? 


0(0) c 


10(3) 


87(3) 


1.61 


1.84 


a 


Possible 7. 1pm band 


4 e 


2.33(3) 


2.74(4) 


0.49(2) 


-0.16(2) 


0(0) 


25(12) 


27(4) 


687(69) 


23(1) 


0.95 


1.15 


b 


Silicate fit too wide 


p Oph E 1 


2.88(3) 


3.36(4) 


0.52(2) 


-0.21(3) 


140(5) 


60(4) 


27(4) 


1206(121) 


98(4) 


0.15 


0.43 


b 


Silicate too wide 


2 


0.23(2) 


0.56(5) 


2.55(12) 


-0.36(5) 


255(16) 


93(7) 


53(8) 


1881(196) 


324(20) 


-0.26 f 


0.51 f 


a 


Deep 7.6/um band 


3 


0.34(38) 


0.65(72) 


1.84(12) 


-0.27(4) 


131(10) 


134(13) 


15(2) 


1536(160) 


155(7) 


-0.06 


0.47 


a 


Very deep 6.8pm band 


4 


0.10(2) 


0.17(4) 


1.42(8) 


-0.23(4) 


291(41) c 


0(0) c 


29(5) c 


1561(279) 


501(49) 


-0.55 f 


-0.01 f 


a 


Very noisy: probably only 
class, COt, r S i good data 


5 


2.62(4) 


5.89(9) 


2.22(8) 


-0.29(4) 


218(14) 


105(7) 


20(3) 


339(34) 


207(8) 


1.17 


1.59 


a 
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4.60(25) 


16.10(89) 


3.31(13) 


-0.52(7) 


306(12) 


117(7) 


0(0) 


193(20) 


112(4) 


1.54 


2.17 


a 




Ser A 1 


0.06(1) 


0.15(3) 


2.14(31) 


-0.60(9) 


845(849)° 


356(161) c 


0(0) c 


2386(311) 


298(33)? 


-0.24 


0.47 


a 


Noisy and faint, possible 
12, 14pm bands 


2 


0.06(1) 


0.17(3) 


2.39(61) 


-0.76(10) 


965(878) c 


298(1 30) c 


0(0) c 


2679(437) 


0(0) 


0.43 


1.35 


a 


Noisy, possible absorption 
bands at 1 1-1 6pm 


3 


0.64(2) 


1.24(4) 


1.70(7) 


-0.32(4) 


523(20) 


165(11) 


19(3) 


1633(163) 


285(13) 


1.35 


1.91 


a 


Possible narrow absorption 
bandat9.7yum-CH 3 OH? 


4 


0.07(1) 


0.15(2) 


2.14(34) 


-0.46(7) 


847(688) c 


271(69) c 


0(0) c 


2426(552) 


269(51)? 


0.42 


1.17 


a 


Noisy, possible absorption 
bands at 11-1 6pm 


5 


0.17(1) 


0.30(2) 


1.42(12) 


-0.37(5) 


435(26) 


165(11) 


8(1) 


1538(177) 


186(11) 


0.25 


0.71 


a 


Possible 12, 14^/m bands 


6 


0.46(1) 


0.45(1) 


-0.18(1) 


-0.22(3) 


272(17) 


97(6) 


0(0) 


803(80) 


223(9) 


0.44 


0.50 


b 


Strong 1 1pm emission 


7 










954(3 12) d 








504(59) 






a 





a Values of r si between -0.8 and 0.6 are considered to be questionable, due to the presence of significant silicate emission and absorption (see discussion in Sect. l3.2.5l . 

b Classification scheme proposed in Sect. l3.2.5l 

c Short wavelength features fitted using linear continuum, fitted from 5. 6-7. 9pm, in cases where the automatic fitting procedure produced errors (see Sect. l2.5> . 

d Feature fitted manually as spectrum was incomplete. 

c Source fell on bad column on array: data from these sources is not considered to be especially accurate. 

f No K-band data: a evaluated using A[ = 8.0pm. 





8yum 


8yum 








Equivalent Line Widths (nm) 




Source 


Cont. 






Source 


flux 

(Jy) 


continuum 
flux (Jy) 




Tllfim 


6/im 
H 2 


6.8/um 


7.6/um 
CH 4 ?? 


13yum 
H 2 


15.2/mi 
CO, 


a 


a 


Group b 


Notes 


SerB 1 


0.19(1) 


0.57(3) 


2.76(16) 


-0.39(5) 


703(44) 


215(31) 


0(0) 


962(99) 


267(14) 


2.02 


2.63 


a 


Very deep absorption features 


2 


0.19(1) 


0.16(1) 


-0.41(2) 


-0.23(3) 


91(5) c 


39(3) c 


25(4) c 


0(0) 


67(10) 


0.14 


0.04 


b 


Large discontinuity at join, 
poor fit from 8-12/jm 


3 


1.37(2) 


0.77(2) 


-1.51(5) 


-0.10(1) 


208(1 6) c 


33(2) c 


43(6) c 


1198(121) 


202(8) 


1.14 


1.05 


c 


Strong silicate emission, 
fit too narrow 


4 


0.33(1) 


0.28(1) 


-0.37(1) 


-0.19(3) 


154(9) c 


45(4) c 


32(5) c 


0(0) 


102(5) 


-1.03* 


-1.14* 


b 


Poss. 13yum absorption band 


5 


0.10(1) 


0.08(1) 


-0.57(3) 


-0.22(3) 


208(8) 


57(4) 


29(4) 


0(0) 


59(4) 


0.06 


-0.09 


b 


Long A fit probably wrong 


6 


6.16(14) 


5.31(13) 


-0.33(1) 


-0.12(2) 


137(5) c 


38(2) c 


24(4) c 


0(0) 


103(4) 


-0.11 


-0.18 


b 


Poor 8-ll//m fit 


7 


0.35(1) 


0.32(1) 


-0.14(1) 


-0.15(2) 


152(7) c 


50(4) c 


23(3) c 


0(0) 


100(4) 


-1.13* 


-1.10* 


b 


Large discontinuity at join 


8 e 


0.31(1) 


0.27(1) 


-0.36(1) 


-0.11(1) 


62(4) 


0(0) 


KD 


554(56) 


27(3) 


-0.31 


-0.29 


b 
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0.35(1) 


0.53(2) 


1.02(4) 


-0.20(3) 


237(8) 


38(4) 


25(4) 


1119(114) 


144(6) 


0.99 


1.33 


a 


Dubious 6.8yum and 7.5/mi fits 


10 


0.36(1) 


0.33(1) 


0.04(1) 


-0.06(1) 


168(1 l) c 


35(3) c 


23(3) c 


609(61) 


126(6) 


-0.43 


-0.34 


b 


Probably peak in nebulosity 


11 


0.36(2) 


1.89(10) 


4.34(31) 


-0.62(8) 


984(45) 


350(67) 


73(11) 


1996(202) 


479(28) 


2.07 


3.14 


a 


Very deep absorption features 


12° 










414(76) c ? 


0(0) c ? 


19(5) c ? 


- 


143(16)? 








Very dubious measurements: 
peaks on bad column 


13 


0.04(1) 


0.07(2) 


1.35(22) 


-0.20(3) 


348(5 l) c 


146(46) c 


69(1 l) c 


2496(291) 


185(7) 


0.55 


1.09 


a 


Good fit but very faint 


Cha I l e 


0.71(2) 


0.65(2) 


-0.18(1) 


-0.07(1) 


0(0) 


0(0) 


7(1) 


548(55) 


0(0) 


-0.32 


-0.28 


b 


Poor fit, due to large 
discontinuity at join 


2 


- 


- 


- 


- 


- 








150(6) 


1.81 


- 


a 




3 


0.14(1) 


0.13(1) 


-0.37(2) 


-0.07(1) 


144(10) c 


36(3) c 


10(2) c 


0(0) 


0(0) 


-0.40 


-0.52 


b 


Noisy, poor fit at long A 


4 


0.04(1) 


0.02(1) 


-1.88(11) 


0.13(2) 


225(25) c 


87(9) c 


0(0) c 


0(0) 


0(0) 


-0.65 


-1.09 


c 


Very faint, low S/N 


5 


3.58(6) 


3.59(6) 


0.00(1) 


-0.12(2) 


41(2) c 


0(0) c 


7(D C 


283(28) 


0(0) 


-0.43 


-0.37 


b 


Strong 1 1/jm emission 


6 


3.18(6) 


1.93(4) 


-1.32(5) 


-0.04(1) 


0(0) 


0(0) 


0(0) 


0(0) 


0(0) 


-0.85 


-1.07 


c 


Silicate fit too narrow 



" Values of Ts, between -0.8 and 0.6 are considered to be questionable, due to the presence of significant silicate emission and absorption (see discussion in Sect. 13.2.51 . 

b Classification scheme proposed in Sect. l3.2.5l 

c Short wavelength features fitted using linear continuum, fitted from 5.6-7.9/mi, in cases where the automatic fitting procedure produced errors (see Sect. l2.5t . 

d Feature fitted manually as spectrum was incomplete. 

c Source fell on bad column on array: data from these sources is not considered to be especially accurate. 

* No K-band data: a evaluated using A[ = 8.0/jm. 
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is broadly consistent with the value of 5.1xl0 17 cm 2 measured 
here (assuming a band strength of A=l.lxl0~ 17 cm molecule -1 , 



Gerakines et al. [1995). Keane et al. (2001 1 measure optical 
depths of T(, .o=0.23 and T6.s=0.09, which are consistent with 
our results of 0.23+0.01 and 0.13+0.01 respectively. It should 
be noted, however, that our value of the 8//m flux density, 
1.05±0.02xlO~ 16 W cm" 2 //m _1 , is somewhat less than the 
value of ~1.7xlO~ 16 W cirT 2 yUirT 1 measured by Whittet et 
al. 0996 1. This may indicate a possible error in the absolute 
flux calibration of the CVF data. However, the measurement of 
spectral features depends only on the relative flux calibration 
from pixel to pixel, which is considered to be accurate through- 
out. 



5 8 10 12 14 15 

A / /Am 

Fig. 4. The spectrum of HH100-IR, with the fitted spectrum 
(solid line) and continuum (dashed line). Note the poor fit to 
the CO2 feature at 15/im as a further example of why the CO2 
feature was fitted independently. 

detected, this was found to be the most reliable and consistent 
method of spectral fitting. 

In order to check the instrumental calibration and the valid- 
ity of this fitting method, comparisons were made with previ- 
ous observations of the bright, well-studied object HH100-IR 
(RCrA 5, see Fig©. Whittet et al. J1996i find a silicate opti- 
cal depth of 1.21+0.05, which is somewhat less than the value 
of 1.35+0.05 measured here: this is almost certainly due to the 
fact that the iterative fitting procedure used here allows the con- 
tinuum to drift away from the observed spectrum slightly. As 
discussed in Sect l2.4l this is due to absorption across the entire 
wavelength range and is a benefit made possible by the broad 
wavelength coverage of the CVF. Whittet et al. ( 1996 1 also find 
a H2O ice column density of 2.4xl0 18 cirT 2 , based on observa- 
tions of the 3/im stretching mode. Adopting a band strength of 
A=1.2xl0~ 17 cm molecule -1 (Gerakines et al. 1995 1, we evalu- 
ate the column density TV as: 

f r(A)dA 



where J T(A)dA is the equivalent width. This gives a value of 
4.8xl0 18 crrT 2 , approximately double the value obtained from 
the stretching mode. This is consistent with the "6/zm to 3/im 
paradox" (Gibb et al. 2000 Dartois & d'Hendecourt l2o6H . 
which commonly results in observations of these two bands 
producing column densities which differ by a factor of ~2. The 
recent work of Gibb & Whittet (2002) suggests that the ex- 
cess depth of the 6/im feature is the result of blending with 
a feature due to organic refractory matter. They also find a 
strong correlation between the excess absorption in the 6/im 
water ice feature and that of the 4.62jum "XCN" feature; un- 
fortunately the wavelength range of the CVF prevented ob- 
servation of the "XCN" feature here. Similarly, from observa- 
tions of the 4.27/vm stretching mode Nummelin et al. l|2001) 
infer a CO2 ice column density of 6.2±0.6xl0 17 cirT 2 , which 



3.2.1. Water Ice: Bending and Libration Modes 

The observed equivalent widths of the 6/mi bending and 13/mi 
libration modes of water ice showed no significant correlation 
at all. Further, the libration mode was not observed to correlate 
with any of the measured features. Whilst there is evidence for 
this in the literature (eg. Bowey at al. 1998 1 the most likely ex- 
planation is that the libration mode is poorly fitted, due to the 
wide variation in possible profiles. The peak of this band has 
previously been found at ~1 ljt/m in crystalline water, ~12.5//m 
in amorphous water and at even longer wavelengths in mix- 
tures with other molecules (Hagen et al. [T9E3i d'Hendecourt 
& Allamandola 1986 Cox 1989). Here a single profile was 
used (that of amorphous ice) and this probably resulted in poor 
fitting. Unfortunately this section of the spectrum is strongly 
blended with the silicate feature, so this problem cannot be 
remedied without making further ad hoc assumptions about the 
silicate profile(s). 



3.2.2. The Unidentified Feature at 6.8//m 

As seen in Fig|5| the strength of the unidentified 6.8/im feature 
is observed to correlate far more strongly with the neighbouring 
6/vm band of water ice than with the 15.2//m band of CO2 ice. 
This is consistent with previous observations, where the 6 and 
6.8yimi features usually correlate well (see Schutte et al. 1996 
and the review by Schutte 1997). This suggests that the carrier 
of the unidentified feature is probably a strongly polar ice, as 
its presence matches far more closely the strongly polar H2O 
ice than the markedly less polar CO2. The more recent work of 
Keane et al. (2001) found the peak position of this feature at 
different wavelengths towards different sources, and proposed 
that it consists of two (related) components. Further detailed 
study of the 6.8yum profiles could yield a great deal more in- 
formation about their chemistry. However, at the low spectral 
resolution of the CVF the 6.8yum profiles are essentially con- 
stant across all the sources, and so such an investigation was 
not possible here. One individual source is also worthy of note: 
GY262 (p Oph E 3) shows a deep 6.8/mi feature with no corre- 
sponding 6.0/im feature (the "rogue" point to the upper left in 
FigHJ. 
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Fig. 5. The strength of the unidentified feature plotted against 
the bending modes of both water (top) and CO2 (bottom) ices. 



3.2.3. 11.2/zm Feature 

The depth of the measured 11.2/xm feature was found to cor- 
relate strongly with the depth of the silicate feature, with a di- 
rect (negative) proportionality providing a good fit to the data. 
Consequently the 1 1 .2/mi feature was interpreted to be an emis- 
sive shoulder on the silicate feature, rather than an independent 
feature due to another species: the presence of this feature nar- 
rows the silicate absorption profile slightly. The silicate profile 
has previously been found to vary depending on the composi- 
tion and structure of the silicate grains (eg. Demyk et al. 2000 1, 
so such a shoulder is not unexpected. However, 3 of the sources 
(RCrA 1, Ser A 6 and Cha I 3) show a significant 1 1.2/rni emis- 
sion feature without any significant silicate absorption feature: 
these features may be attributable to emission from crystalline 
silicates (Bregman et al. 1987 , Campins & Ryan 1989 1. 

3.2.4. C0 2 Ice Profiles 



As mentioned in Sect. 12.41 the CO2 ice profile is a very sensi- 
tive diagnostic of the ice environment. Whilst the spectral res- 
olution of the CVF was not sufficient to study the profiles in 
great detail, it is worth noting that a significant long wavelength 
wing was present in almost all the observed absorption profiles. 
This would seem to indicate that a large fraction of the CO2 ice 
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Fig. 6. Examples of observed silicate profiles: emission profiles 
are generally broader than those seen in absorption and com- 
posite emission/absorption profiles are also seen. "Type" refers 
to the classification scheme described in Sect. 13.2.51 

observed exists in a polar (F^O-rich) phase, which is charac- 
terised by this long wavelength wing (Gerakines et al. 1999 1. 
However quantifying the relative abundances of the different 
phases was not possible at this low spectral resolution. 

3.2.5. Spectral Classification and Silicate Profiles 

If we look at the shapes of the spectra, it appears natural to 
divide the sources into 3 distinct groups based on the strengths 
of the spectral features observed. An example of each type is 
shown in Fig[6]and the classifications are included in Table|2] 

a) The first group show deep absorption features due to all of 
the ices and also due to the silicates (typically these have 
ts; > 0.6). These objects are interpreted to be heavily em- 
bedded objects, showing strong absorption due to cold fore- 
ground material. 

b) The second group shows weaker ice absorption features and 
weak silicate features (-0.8 < rs; < 0.6) which are, in 
general, not especially well fitted. These are interpreted as 
less heavily embedded objects, showing less foreground ab- 
sorption. The poor silicate fits are probably due to compli- 
cations caused by combined silicate emission and absorp- 
tion. 

c) The third group show strong silicate emission < -1.0) 
and little or no ice absorption. It should be noted that in 
the 2 good S/N spectra showing strong silicate emission, 
the observed emission profiles were significantly broader 
than both the observed absorption profiles and the profile 
used to fit the data (see Fig|6). (However it should also be 
noted that one of these 2 sources, SVS2, is an unresolved 
double.) These are probably young stars which have shed 
most of their circumstellar matter, retaining only a little hot, 
optically thin silicate dust around them (possibly in a disc). 
It therefore seems plausible that these three groups might 
represent an approximate evolutionary sequence, with the 
YSOs shedding circumstellar material as they evolve. 
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Fig. 7. Composite silicate profiles fitted to the spectrum of 
SVS20. The solid line shows a composite fit with r a bs=0.59 
and T em =-0.33, the dashed line a composite fit with T a bs=0.88 
and T em =-0.48 (and a different continuum). The single profile 
fit reported in the results table, with ts;=-0.33, is shown as the 
dotted line. 



In order to test this the spectrum of SVS20 (Ser B 6) was 
re-fitted using a composite silicate profile, consisting of the 
standard absorption profile used above superimposed on the 
broader emission profile of S VS2 (Ser B 3 - see Fig®. As seen 
in FigQthis provided a good fit to the observed silicate profile, 
indicating that both silicate emission and absorption probably 
do occur, with different profiles seen in absorption and emis- 
sion. The physical interpretation of this is that these objects 
have a hot region of optically thin silicates around them but 
remain embedded in the cloud, which results in foreground ab- 
sorption due to silicates and ices being superimposed on the 
"intrinsic" silicate emission feature in the spectra. 

It should be noted that while it is possible to fit such profiles 
to these sources, they are not well constrained. As can be seen 
by comparing the 2 composite profiles in Fig0there is a degen- 
eracy between the strength of this "flat-topped" silicate profile 
and the continuum strength, and a great deal of information is 
assumed about the silicate profiles themselves. Consequently 
this investigation was not pursued for all the sources. Where 
-0.8 < T$i < 0.6 the values in Table [2 represent the correct 
optical depth at 9.7/im, but are usually in error towards the 
wings of the silicate feature. Such values are useful in studies 
of general trends, but they do not accurately describe the en- 
tire silicate feature: a single number cannot fully describe the 
complicated nature of these features. 

This also provides more direct evidence, similar to that of 
Whittet et al. ( 1988 ), that the extensive scatter in measurements 
of the Ay/rsi ratio towards YSOs may be caused by complica- 
tions due to silicate emission. Direct measurements of the sili- 
cate optical depth will significantly under-estimate the depth of 
foreground material in cases (such as SVS20) where optically 
thin silicate emission is also present. Rieke & Lebofsky ( 1985) 
derive Av/ts; - 17, whilst noting the problem posed by "in- 



trinsic" silicate emission, and such emission does indeed lead 
to errors when evaluating Ay from the silicate depth. As an ex- 
ample, the ISOCAM survey of Bontemps et al. (2001 ) derives 
values of A v * 28 for GY252 (p Oph E 1) and * 24 for GY262 
(p Oph E 3); the X-ray survey of Imanishi et al. (2001 ) derives 
slightly lower values of ~ 24 and ~ 1 8 respectively. However 
GY252, as seen in Fig|6] clearly shows a "composite" silicate 
profile whereas GY262 does not, and the respective silicate op- 
tical depths of 0.52 and 1.84 are clearly inconsistent with a 
single Ay/rsi ratio. Unless the emitting and absorbing compo- 
nents of such silicate features can be separated unambiguously 
obtaining a single Ay/rs; ratio will not be possible. 

3.2.6. Spectral Index 

Of the 41 objects to which spectral indices were assigned 20 
were found to be class I objects (ffcont > 0) and 17 to be class 
II (-1.5 < a cont < 0). In addition to this 3 objects (p Oph E 
1 and 3, Ser B 2) were found to lie very close to the class I/II 
border (so-called transition objects: Andre & Montmerle l 1 9941 
Greene et al. 1994). The lack of any class III objects is a se- 
lection effect, but it seems to indicate that all of the observed 
sources are within the clouds, as field stars would show little 
or no infrared excess and would appear as class III objects in 
this survey. The distribution of class I and II objects is approx- 
imately constant throughout all four regions, and no significant 
trends involving spectral class were found. Evaluation of mid- 
IR spectral indices appears to support the finding of Bontemps 
et al. (2001) that mid-IR spectral indices do not discriminate 
strongly between class I and class II objects, and a classifica- 
tion scheme similar to that in Sect. l3.2.5l seems preferable when 
considering spectroscopic data. 

Comparing the index-based spectral classes to the spectral 
feature grouping above, we see that the majority of the class 
I sources belong to group a) and the majority of the class II 
sources to group b). However it is notable that a significant 
number (6) of the 23 group a) objects are either class II or tran- 
sition objects. This in keeping with the convention that class I 
objects are the youngest (Lada 1987 1, but it may be significant 
that there are two distinct types of spectrum which both result 
in class II SEDs. Objects showing strong foreground absorp- 
tion and objects showing little foreground absorption can both 
present class II continuum SEDs. This could be due to geom- 
etry of the YSOs relative to the absorbing clouds, or could be 
an effect intrinsic to the sources: spectroscopy along a pencil- 
beam cannot distinguish these. Whilst a larger sample size 
would constrain this problem further, it does seem that mid- 
IR observations alone do not provide strong constraints on the 
"true" spectral class of YSOs. 

4. Variations Between Target Regions 

In general, the most striking regional variation is between the 
sources in Cha I and those in the other three regions. All of the 
sources in Cha I, except for ISO-Chal 192 (Cha 1 2), show little 
or no absorption due to volatile ices and show silicates either 
in emission or an emission/absorption composite. As it lay out- 
side the CVF1 images, due to rotation of the spacecraft between 
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Fig. 8. The strength of the bending mode of water ice plotted 
against the silicate optical depth. Different regions are distin- 
guished by different symbols, as shown in the legend. 



Fig. 9. The ratio of the CC^I^O ice equivalent widths plotted 
against silicate optical depth. This ratio is directly proportional 
to the ratio of ice column densities: adopting the band strengths 
from Gerakines et al. ( 1995 1 results in an equivalent width ratio 
of 1 corresponding to a column density ratio of ^0. 17. 



the two sets of observations, we only have a 9.3-16.3/im spec- 
trum for ISO-Chal 192. However it shows CO2 ice absorption 
and also appears to show deep silicate absorption (rs; ~ 2—4), 
so it is probably an embedded object. Persi et al. ( 1999) found 
that ISO-Chal 192 lies in a dense core and that it is probably the 
engine for the bipolar outflow detected by Mattila et al. £19891. 
Jones at al. £1985) found that several other YSOs have formed 
around the edge of this core, rather than at its centre: these 
were apparently formed during a burst of star formation trig- 
gered by a wind from the nearby star HD97300. Our results are 
consistent with this: ISO-Chal 192 appears to be heavily em- 
bedded, and has a very red SED. The rest of the surrounding 
objects show similar spectral characteristics, and appear to be 
surrounded by far less circumstellar/foreground material than 
ISO-Chal 192. 

12 of the 20 objects in Serpens show deep absorption 
features, 7 show composite silicate profiles, and one (SVS2) 
shows silicate emission. It is clear that Serpens contains a wide 
range of YSOs, ranging from the most heavily embedded ob- 
ject observed (Ser B 1 1) to a source showing silicate emission. 
By contrast almost all of the sources in both RCrA and p Oph 
show deep absorption features, with only 1 of the 7 sources in 
RCrA and only 2 of the 10 sources in p Oph showing compos- 
ite silicate profiles. Given that YSOs are generally expected to 
"sweep out" their circumstellar material as they evolve (Shu et 
al.[T987i, this seems to imply that the observed star formation 
is at a similar evolutionary stage in both RCrA and p Oph, prob- 
ably an earlier stage than that in Serpens. However, given the 
broad variation in the sources observed in Serpens, this conclu- 
sion remains somewhat tentative. It should also be noted that 
the observed fields (of 2'x2') are much smaller than the star- 
forming clouds, which typically extend over several square de- 
grees of the sky, so the populations observed here may not be 
representative of the entire clouds. 



4.1. Silicate Depth and H 2 Bending Mode 

Figure[S]shows the measured equivalent widths of the 6pm wa- 
ter ice feature plotted against the optical depth of the silicate 
feature: clear regional variations are apparent. As discussed 
above, all but one of the sources in Cha I are characterised by 
little or no ice absorption and by silicate emission; the other 3 
regions show absorption from both silicates and ices with only 
a few examples of weak silicate emission. The sources in RCrA 
are characterised by similar absorption strengths throughout, 
whereas far deeper absorption and variation is observed in both 
p Oph and Serpens. It is clear from Fig|8] however, that the 
sources in p Oph show markedly less water ice absorption rel- 
ative to silicates than those in Serpens. This could be due to 
sublimation effects, indicating a difference in the thermal struc- 
ture and history of the two clouds, or could represent a real 
difference in the chemical composition of the two clouds; with- 
out measuring the corresponding gas-phase lines we cannot say 
which. Sub-millimetre observations of gas-phase H2O (Ashby 
et al. 2000 1 do show an under-abundance of gaseous H2O in 
p Oph A relative to S140, but to the best of our knowledge 
no similar direct comparisons between the sources observed 
here exist in the literature. It is also worth noting that the sili- 
cate:water ice ratios appear to be approximately constant for all 
the different embedded sources in both p Oph and Serpens. This 
indicates that in both clouds the composition of the absorb- 
ing intracloud material is roughly constant along different lines 
of sight, a result consistent with the previous Serpens work of 
Eiroa & Hodapp ( fT9"S9l . 

4.2. H 2 and C0 2 Bending Modes 

Figure|9]shows the ratio of the C02:H20 ice equivalent widths 
(which is proportional to the corresponding ratio of column 
densities) plotted against the silicate optical depth. The re- 
sults from Serpens and RCrA are fairly well correlated, with 
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CC>2:H20 ice column density ratios ranging from to 0.16 
(adopting the band strengths from Gerakines et al. 1995 1. There 
is a possible trend showing the CCVF^O increasing with de- 
creasing silicate optical depth, but this is somewhat unclear in 
these data. 

The data from p Oph, however, show a far greater scat- 
ter with CC>2:H20 ice column density ratios ranging from to 
0.4. This is consistent with previous observations: for example 
Chiar et al. i 19941 [1995 ) found a greater scatter in column den- 
sities of CO and H2O ices against A y in p Oph than in Serpens, 
RCrA or Taurus. The reasons for this are unclear but it seems 
likely that, as suggested by Chiar et al. J1995I I. local conditions 
around the young stars in p Oph play a role. Detailed study of 
the gas-phase abundances of CO2 and H2O towards these ob- 
jects will further constrain this problem, but without such ob- 
servations or further knowledge of the envelope structures we 
merely note the presence of this discrepancy here. 

5. Conclusions 

ISO-CVF spectroscopy of four low-mass star-formation re- 
gions has been performed, and a total of 42 different sources 
were observed. The resulting low-resolution 5-16. 5/im spec- 
tra have been analysed, and a number of different conclusions 
drawn: 

- The profile of the 9.1pm silicate absorption feature varies 
somewhat from source to source, but usually shows a weak 
emissive shoulder at ^ 1 1 .2pm. Significantly broader pro- 
files are observed in emission. A number of sources appear 
to show "composite" profiles, with absorbing and emitting 
components superimposed. 

- The presence of a long wavelength wing on the 15.2jum 
CO2 ice band in most of the sources suggests that a sig- 
nificant fraction of the CO2 ice observed exists in a polar 
(H20-rich) phase. 

- The strength of the unidentified 6.8//m band is observed to 
correlate far more strongly with the 6.0pm water ice band 
than with the 15.2yum CO2 ice band. This suggests, in a 
manner consistent with previous observations, that the car- 
rier of this band is a strongly polar ice. 

- Compared to the other 3 regions (which show similar spec- 
tral characteristics) the sources observed in Cha I are some- 
what anomalous, showing comparatively little foreground 
absorption. However this is consistent with previous stud- 
ies, which indicate that the observed Cha I sources repre- 
sent an unusual cluster of YSOs which have formed at the 
edge of a dense cloud core. 

- The p Oph cloud appears to have a systematic under- 
abundance of water ice, relative to silicates, relative to both 
RCrA and Serpens. 

- The C02:H20 ice ratios towards the sources in p Oph show 
considerably greater scatter than in the other regions, with 
values ranging from to 0.4. This may be due to local con- 
ditions around the YSOs in p Oph but more data is needed 
to constrain this problem further. 

The most important limiting factor in this study is obviously 
the low spectral resolution of the CVF, which essentially pre- 



vents any detailed study of absorption profiles, particularly 
around the CO2 ice feature. However the use of the CVF en- 
abled the study of a large number of objects, many of which 
would not have been observable at higher spectral resolution 
due to the low S/N. Consequently we have obtained a unique 
data set, with a large sample size and broad wavelength cov- 
erage. This enabled simultaneous study of the dominant mid- 
IR spectral features along the line-of-sight to a large number 
of YSOs, in a manner not possible using other spectroscopic 
methods. Further detailed study of the ice absorption profiles 
will undoubtedly yield more information regarding the chem- 
istry around these objects, and study of the silicate profiles 
may yield more information about the dust composition and 
the structure of the dust envelopes. 
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Appendix A: Atlas of Spectra 

All of the spectra are presented in this appendix. In each plot 
the data points are shown as crosses, with the fits as solid lines 
and the fitted continua as broken lines. The spectra are arranged 
in the same order as in Tabled The error bars shown are the sta- 
tistical errors on the data points and do not include the absolute 
flux calibration errors discussed in Sect. 13.21 Please also note 
that the fits shown are the "complete" fits, covering the com- 
plete wavelength range: in some cases (as indicated in Table 
13 and discussed in Sect. 12. 5> the narrower features were fitted 
independently and these fits are not shown here. 



